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Abstract: 



We propose that in QCD with dynamical quarks, colour deconfinement occurs when 
an external field induced by the chiral condensate strongly aligns the Polyakov loop. This 
^ ' effect sets in at the chiral symmetry restoration temperature and thus makes decon- 

^ ' finement and chiral symmetry restoration coincide. The predicted singular behaviour of 

^ . Polyakov loop susceptibilities at is shown to be supported by finite temperature lattice 

! calculations. 

O 

o 
o 

QhI Finite temperature QCD at vanishing overall baryon density leads to two well-defined 

Q^l phase transitions. In the limit of infinite bare quark mass, for rriq oo, one obtains pure 

SU{N) gauge theory; here deconfinement occurs when the global center G SU{N) 
symmetry of the Lagrangian is spontaneously broken. The expectation value L{T) of 
the Polyakov loop constitutes the order parameter for this transition, analogous to the 
^ I magnetization m(T) in Z^r spin theories. There exists a critical temperature T^, with 

d I L{T) = V T < Tfi and L{T) > V T > T^. The introduction of dynamical quarks 

{rriq < oo) breaks this symmetry explicitly, and now L{T) > V T > 0. 

For nig = 0, the chiral symmetry of the Lagrangian is spontaneously broken at low 
temperatures and restored for T > T^. The chiral condensate K{T) = (ipip) is the order 
parameter for this transition, with K{T) > V T < and K{T) = V T > T^. 
A non- vanishing quark mass breaks chiral symmetry explicitly, and for large rrig, the 
temperature variation of the chiral condensate becomes completely smooth. 

One might consider the inverse bare quark mass to play such a symmetry breaking 
role for deconfinement, which would imply that for — 0, the temperature variation 
of the Polyakov loop expectation value would also become smooth In lattice studies 
it is found, however, that this is not the case p|: even for uig 0, the Polyakov loop 
varies sharply with temperature and the corresponding susceptibility {{L"^) — (L)"^) peaks 
sharply at the chiral restoration temperature T^, which is considerably lower than T^. 
In somewhat loose terminology, one describes this situation by saying that in QCD for 
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niq — > 0, deconfinement and chiral symmetry restoration coincide. The aim of this paper 
is to elucidate the underlying reasons for such behaviour and to show that chiral symmetry 
restoration in fact drives the observed Polyakov loop variation; in particular, it leads to 
singular behaviour for specific Polyakov susceptibilities. 

Conceptually, our starting point is the idea that it is the inverse of the constituent 
quark mass tuq, rather than the bare mass m^, which acts as an external field for the 
Z]\f symmetry [0]. Evidently there does not exist a clear definition of mg; for the present 
discussion, we take it to be determined by the (non-Goldstone) hadron masses. For our 
actual conclusions, however, we shall avoid the intermediate step of constituent quarks 
altogether. In the limit rriq —>■ 0, itlq remains finite in the temperature region in which 
chiral symmetry is spontaneously broken, since the hadron masses do. If the hadron 
masses are related to the chiral condensate, e.g., if the nucleon mass is given by [Q 

M„ ~ 3 mg ~ (^^)^/^ (1) 

then chiral symmetry restoration with (tptp) —>■ will lead to a sudden increase of the 
external field i7 ~ 1/ mq, forcing a large explicit breaking of Zn and hence a corresponding 
alignment of the Polyakov loop. We want to argue that it is this effect which causes 
a sharp variation of the Polyakov loop at < T^, with the functional form of the 
variation determined by chiral symmetry restoration, which is in general different from 
that obtained in pure gauge theory at T^. 

We begin by summarizing the main features of the critical behaviour for the two 
transitions, based on QCD with two species of dynamical quarks and colour SU{3). The 
order parameter for chiral symmetry restoration, K = vanishes for niq = as 

T ^ from below, 

if(T,m, = 0)~(T^-Tf, T<T^ (2) 
while at T = T^, it vanishes as 

K{T^ = 0,mq)r^ml/' (3) 
when rriq 0. The chiral quark mass susceptibility 

x^(T,m,) = (^^^^~ ~ \T-T,r (4) 

diverges at T = for rrig —>■ 0. The critical behaviour at the chiral symmetry restoration 
point is thus specified in terms of critical exponents j3, 6, 7, .... For two massless quark 
flavours, the transition is conjectured to belong to the 0(4) universality class 0, which 
would determine the exponents. Present lattice calculations cannot yet determine if this 
conjecture is correct 0, |^. 

In the limit rUq — > 00, we recover pure SU{3) gauge theory, for which the deconfinement 
transition is of first order, as it is for the corresponding Z3 spin theory 0. The resulting 
discontinuity at T = persists in L(T, nig) for sufficiently large m^; but below some quark 
mass value < < 00, L(T, m^) and all its derivatives presumably vary smoothly with 
temperature. We want to argue, however, that for rriq 0, certain derivatives of L{T, rriq) 
will become singular. 
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As indicated above, we now take the Z^r structure of QCD to be that of Z^v spin theory 
in an external field if, determined by the constituent quark mass. From Eq. (|l|), H and 
hence also the Polyakov loop thus become functions of K: the quark mass dependence 
of L enters through K{T,mq), with L{T,K). In spin systems with external field, the 
magnetization m(T, H) is for H ^ an analytic function of T and H. We therefore 
assume that L{T, K) is for oo > > an analytic function of T and K, with the quark 
mass dependence of L entering through K{T,mq). 

As a consequence, the critical behaviour of the chiral condensate will be reflected in 
the behaviour of the Polyakov loop at T = T^. A first hint that this is the case is seen in 
|P-|Tl|. For four flavours of light quarks, the chiral transition becomes first order, with 
a discontinuity in K{T,m) at for ruq < smaller than some 'endpoint value' m^. 
The corresponding Polyakov loop L{T) also shows a discontinuity at this temperature, 
as expected from our considerations. Note that here we have a first order transition at 
T = for < rUq < m^, induced by chiral symmetry restoration, and another first order 
transition at T = Td > for rriq > rriq, induced by spontaneous Zn symmetry breaking. 
For < rriq < rriq, L{T, K) varies smoothly with T. 

In the case of two dynamical quark flavours, the chiral transition is presumably con- 
tinuous. From the L-variation 

we then find that the Polyakov loop susceptibilities 

(6) 





drriq y y Y dK J ^\ drriq j ^ 

and 




must diverge at T = T^ in the chiral limit = 0, since the chiral susceptibilities Xm — 
{dK/dmq)T and Xt = {dK/dT)mg diverge in this limit. 

In 0, the chiral susceptibilities Xm Xt were studied on an 8'^ x 4 lattice for 
quark masses rriqa = 0.075, 0.0375 and 0.02. In Figs. |1] and |^, the results are shown as 
functions of the effective temperature variable k = 6/(7^, where g denotes the coupling in 
the QCD Lagrangian. The increase of the peak height for decreasing quark mass indicates 
the divergence in the chiral limit rriq — 0. In Figs. |^ and |, we show the corresponding 
results for the Polyakov loop susceptibilities. They also peak sharply, the peak positions 
coincide with those for the chiral susceptibilities, and here as well the peak height increases 
with decreasing quark mass. The observed behaviour therefore provides clear support for 
the divergence of the temperature and quark mass derivatives of the Polyakov loop at 
in the chiral limit. 

A direct comparison of the functional behaviour, ideally of the relevant critical expo- 
nents, for the two cases becomes difficult for two reasons. The chiral susceptibilities do not 
at present lead to the predicted 0(4) exponents; this may indicate that the quark masses 
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are still too large for pure critical behaviour. The Polyakov loop susceptibilities contain in 
addition unknown non-singular factors {dL/dT)^ and {dL/dK)T, which will modify the 
non-singular m^-dependence relative to that of the chiral susceptibilities. Lattice studies 
for smaller quark masses would therefore be very helpful. 

The 'temperature' susceptibilities xf and Xk have also been studied on larger lattices 
12^ X 4 and 16^ x 4, for the same three quark mass values [§. In Figs. |^ and ^ we show 
the dependence of the peak height on the spatial volume. The two susceptibilities show 
a very similar peak increase in going from 8'^ to 12'^, but then saturation, indicating that 
there are no further finite size effects for the given quark mass value. 

In summary, we have shown that chiral symmetry restoration, with the resulting 
sudden change in an effective constituent quark mass, leads to a suddenly increasing 
external field which aligns the Polyakov loops at T = < and thus produces a strong 
explicit breaking of the Zn symmetry of the gauge field part of the Lagrangian. This 
effect makes chiral symmetry restoration 'coincide' with deconfinement. The resulting 
predictions for diverging Polyakov loop susceptibilities are found to be well supported by 
finite temperature lattice calculations for full QCD with two flavours of light quarks. 
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Figure 1: The chiral temperature susceptibility Xk function of the temperature variable 



30.0 



20.0 



10.0 



0.0 




Figure 2: The chiral quark mass susceptibility Xm function of the temperature variable 




Figure 3: The Polyakov loop temperature susceptibility Xk function of the temperature 
variable k = Q/g^. 




mq=0.02 



5.25 5.30 5.35 5.40 

K 



Figure 4: The Polyakov loop quark mass susceptibility Xm function of the temperature 
variable n = Q/g"^. 
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